Optical supercavitation in soft-matter 
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We investigate theoretically, numerically and experimentally nonlinear optical waves in an absorb- 
ing out-of-equilibrium colloidal material at the gelification transition. At sufficiently high optical 
intensity, absorption is frustrated and light propagates into the medium. The process is mediated 
by the formation of a matter-shock wave due to optically induced thermodiffusion, and largely 
resembles the mechanism of hydrodynamical supercavitation, as it is accompanied by a dynamic 
phase-transition region between the beam and the absorbing material. 



Nonlinear optical propagation in complex fluids is me- 
diated by several and often competing effects, which are 
in many respect un-explored, specifically in the pres- 
ence of structural phase transitions. Among the var- 
ious involved phenomena, there are light induced re- 
orientational effects [lHH , electrostriction_JlJ3 , thermal 
[§-[l2j and thermodiffusive phenomena Hf 15j . 

Theoretical and experimental works have recently put 
emphasis on electrostriction and the corresponding soli- 
tary waves and beam self-trapping [l6l - [22j |. No evidence 
is reported on optical nonlinear waves in the presence of 
typical complex processes, such as dynamic phase tran- 
sitions and aging (see, e.g, 23-2^]). Self-induced trans- 
parency with respect to scattering losses in a non ab- 
sorbing medium, in absence of structural transitions and 
thermal effects, was predicted in ref. 271 ]. 

In this work we consider the nonlinear optical propa- 
gation of a micron-sized focused beam in a dye-doped, 
strongly absorbing, out-of-equilibrium nanoscale colloid, 
which undergoes a dynamical arrest, i.e., a liquid-gel 
transition. We show that, when this happens, the sys- 
tem behaviour changes from being strongly absorbitivc, 
to a regime that allows beam propagation. This is ac- 
compained by a shock-wave-front in the material that 
forms a depletion region of the absorbing molecules. The 
process resembles cavitation and supercavitation hydro- 
dynamic phenomena (see, e.g., [28]): the shock-depleted 



region has vanishing concentration of colloidal beads and 
constitutes a non-arrested layer. This layer allows light 
propagation through the otherwise gelified and absorb- 
ing medium. This happens in the very same way that a 
shock-induced gaseous region favors propulsion in liquid 
supercavitation. The effect is investigated theoretically, 
numerically and experimentally, and the findings under- 
line the main role of time and s pac e nonlocality [29( in 
enhancing shock- phenomena [ill. l30j] . 

Theory — We consider a colloid able to absorb light 
and undergoing a dynamic transition, from a liquid 
regime to a gelified state (high viscosity). The material 
is formed by a host liquid (water) in which interacting 
nano-sized dye-doped (hence, light absorbing) colloidal 
beads are dispersed. In the presence of strong light ab- 



sorption, the leading phenomena are thermally induced, 
and re-orientational effects and electrostriction are negli- 
gible. Specifically, we consider the variation of the refrac- 
tive index due to the thermal density variation in the host 
liquid, and the Soret effect, which is the diffusion of the 
colloidal beads due to a non-uniform temperature profile. 
As the colloidal particles are driven by the light-induced 
thermal gradient, the local absorption and temperature 
change, so that the real and the imaginary parts of the 
optical susceptibility are both affected. Such a mech- 
anism resembles the photorefractive beam nonlincarity 
[3lj ]. where the optical beam generates a space-charge 
field that, through the electro-optic effect, changes the 
refractive index. In our case, the role of the space-charge 
field is played by the temperature, which drives the par- 
ticle motion. 

Non-equilibrium thermodynamics provide the general 
equations for modelling thermodiffusion. [32j Relevant 
variables are the temperature T and the colloidal con- 
centration c, governed by the coupled fluxes Jt and J c 



d t c = V • J c , d t T = V • J T + D : 



The currents are given by 



J T = D T VT + D D Vc- J c = D c Vc + S T D c c{\ - c)VT, 

(2) 

where Dp is the Dafour coefficient (Dp = hereafter), 
Dt is the thermal diffusion coefficient, D c is the mass 
diffusion, St is the Soret coefficient at the working tem- 
perature To, kx is the thermal conductivity and ao(c) 
the light absorption coefficient detailed below. We ne- 
glect the temperature variation of the Soret coefficient. 
The spatial distribution of the continuous wave optical 
field is given by the paraxial equation 



■ 9A 1 „- , 



k ^<El A = -ia {c)A, (3) 



with \A\ 2 the optical intensity, k — 2irn/X the wavenum- 
ber, n is the refractived index, A the wavelength. An 
is the optically induced refractive index perturbation, 
due to the temperature effect An = ATdn/dT with 
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AT = T — T and ao(c) — «b + cda/dc is the loss co- 
effient, with as the residual host liquid absorption in 
absence of the colloid c = (drn < and d c n > 0). 

Introducing p = c/cq with cq the background con- 
centration, 9 = (T — To) /To the normalized tempera- 
ture, and scaling spatial and temporal variables such 
that £ = x/wo, a = y/wo, C = z/zq, r = t/to, with 
zq = kwQ the diffraction length (wq is the beam waist) 
and to = Wq/Dt, Eqs.([I][3]) are written as 

= V 2 fl + (p + po)^! 2 ; a r p = t/V 2 /? + sV e • (pV e 0)(5) 

with a = ZQCod a n, 6 — k\dTn\ToZo/n, V = Vf !<T + C e ^C' 
e = 1/kzQ = 10~ 2 (C is the unit vector of the £ di- 
rection), and ip — ^4/^4-0 with A\ — kTo/aow 2 ,. In 
©, Po = io a s^o/^co9o,n accounts for the residual 
light absorption in the absence of the colloidal beads. 
i] = D c /Dt measures the strength of particles diffusion 
in terms of heat diffusion (Dt), and decreases (77 << I) 
at the dynamics slowing down, s — StTqDt / D c is the 
normalized Soret coefficient. 
As previously investigated 



14j . the temporal scale for 



the onset of a stationary temperature profile is much 
faster than for thermal diffusion. This is a signature of 
soft-matter, where various independent temporal scales 
intervene. In our case, p is slowly varying with respect 
to 9 and follows adiabatically; analogously ip follow both 
p and 9. 

9 is expressed in terms of p and t/>| 2 , by using Eq.© 
in the stationary regime {d T 9 = 0), and Eq.© becomes 
(po = for simplicity) 



9 P V72 



s(V e p)-(V e 9)-sp 2 W 



(6) 



In the initial stage, p = 1, the temperature gradient term 
V e in ([6]) is the heat flux generated by the absorbed in- 
put beam (centered at £ = and £ = a = 0), which is 
given by 9 = P/^a^J (er) 2 + £ 2 where P = J °° \ip\ 2 rdr 
is the normalized beam-power and r = ^/£ 2 + tr 2 ; as 
obtained by the Green function for the heat equation. 
The beam decays much faster in space than 9 and can 
be treated in §5§ as a Dirac delta (highly nonlocal ap- 
proximation for 9). Correspondingly, in ([6]), the drift 
term is found to scale as P/L 2 where L is the distance 
from the peak input optical intensity. Thus the tem- 
perature gradient induces a drift velocity of the density 
p, such that the density travels faster in proximity of 
the optical field, and slower in the surrounding region. 
To address the resulting formation of a shock, we con- 
sider the one dimensional case (d^ — d a =0), such that 
V 2 6> = e 2 d 2 9 = -p\4>\ 2 and e 2 d c 9 = -f c pl. The heat 
source p\tp\ 2 after the depletion stage is a narrow local- 
ized regions (see FigfJ] below) , where the intensity \ip\ 2 



decays as exp(— 2pa() = exp(— 2aQ, being p = 1. Cor- 
respondingly, e 2 c\# = — pl / (2a) and Io the input peak 
intensity. The resulting equation for p becomes 



sL 



d T p+—pd c p = rje d ( p, 



(7) 



which is the Burgers equation predicting the formation 
of a shock front in the £ direction at the dynamic arrest 
(77 -»• 0) (33}. The term p 2 \ip\ 2 in flBJ} is exponentially 
small due to the light absorption and neglected. Extend- 
ing these arguments to the three-dimensional case leads 
to the prediction of a shock front with an semi-elliptical 
shape in the (r, £) plane, for a sufficiently high power and 
at the dynamic transition for the soft-colloidal matter. 
We stress that for low enough power, and for a matter 
relaxation time comparable with the thermal one (i.e., 
in the liquid phase, r) = 1), the shock is smoothed out 
by the "dissipative" term rje 2 d 2 p in Q. It is this light 
induced shock that produces a depletion of the absorb- 
ing material, which correspondingly becomes transparent 
and allows light propagation. We stress that our model 
describes nonlinear optical propagation both in the liq- 
uid and in the arrested phase, which are described by dif- 
ferent values for r\ (vanishing in the arrested phase and 
determined by the initial concentration cq), i.e., by differ- 
ent values of the matter diffusion coefficient D c , which is 
the parameter mainly affected in the gelification process 
(see also 0,E1)- 

Spatio-temporal dynamics — We numerically solved 
Eqs.(j4][5|) by using a finite-difference time-domain 
predictor-corrector scheme for the temporal evolution of 
p and 9, and a pseudo-spectral beam-propagation scheme 
for i/j. At each instant the distribution of the field is cal- 
culated from the distribution of p and 9. To reduce the 
computational effort, we limited to a one-transverse di- 
mension £ (d a = 0). Parameters are chosen from realis- 
tic values. For an input Gaussian beam (with amplitude 
A and unitary normalized waist), we report representa- 
tive simulations. For a large set of parameter values, 
the overal spatio-temporal dynamics are qualitatively the 
same, as described in what follows. 

After an initial stage for the formation of the tempera- 
ture profile, the thermally driven diffusion starts to act on 
the p field and induces a depletion region (where p « I), 
as in FigfT] As the density reaches the value p = 0, the 
shock-front starts moving in the material. In FigfJ] we 
show three snapshots of the density profile. Figj3] shows 
the longitudinal (£ = 0) and transverse (C = 0) shock 
profiles, superimposed to the temperature, which unveil 
the traveling front driven by 9. Correspondingly, the 
material becomes transparent to the optical field, whose 
propagation is forbidden in the linear regime by strong 
absorption. These dynamics are obtained only at small 
values of 77, while at higher values (77 = 1) the diffusion 
of the material smooths out and rapidly dissipates the 
shock front. FigfJ] shows the spatial beam intensity for 
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FIG. 1: (Color online) Initial stage of the formation of the 
cavitation. Thick line is material density (£ = 0, £ > 0), thin 
line is temperature. Quantities are scaled to their peak value 
(a = 6, 5 = 1, s = 10, ?/ = 0.1, p = 1, A = 0.1). 
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FIG. 2: (Color online) Density distribution p (depletion layer) 
at different instants. Parameters as in Fig[T] 



three different instants. As the matter shock front ad- 
vances, the beam enters the material that becomes pro- 
gressively transparent. We also report the fluorescence 
signal, given by ( o|'0| 2 , which appears as a propagating 
localized wave front and is experimentally accessible. 

Experiments - In our experiments we use clay 
(Laponite) in water, doped by Rhodamine B (RhB), a 
dye-doped colloidal solution displaying a gelification pro- 
cess [3J-|36j. The characterization of the nonlinear sus- 
ceptibility was reported in [l3j, E3) ■ RhB molecules ad- 
sorb on the surface of the Laponite platelets (disks with 
diameter of the order of 25nm), and the light absorption 
depends on c. The sample (length 1cm) is illuminated by 
a focused (wq = 10/Ltm on the input facet) linearly polar- 




FIG. 3: (Color online) Temperature (thin line) driven shock 
front in the longitudinal (£ = 0) and transverse (£ = 0) sec- 
tions (thick line is density). Parameters as in FigfTJ 




FIG. 4: (Color online) Self-transparency mediated by the 
shock front. Top panels: beam intensity. Bottom: fluores- 
cence signal p|"i/>| 2 . Parameters as in Fig[T] 



ized Gaussian beam (frequency doubled continuos wave 
Nd:Yag A =532nm, n = 1.3), obtained through a 40x mi- 
croscope objective. Absorption losses are ao — 5mm -1 
(when c = 2wt%): the sample appears red and is com- 
pletely absorbing the 532nm wavelength. We studied 
two different conditions of concentration cq = 1% and 
co = 2%. For co = 2wt% the parameters previously de- 
fined are z = 1.5mm, dn/dT Kr 4 (°C)-\ T = 300K, 
e 1(T 2 ,(5 10, a = 10 (D C /D T = V = 0.1). The flu- 
orescence signal is imaged by a microscope from the top 
of the sample. 

At low concentration (co = lwt%), the Laponite is 
far from the dynamic slowing down (liquid phase), and 
even at very high power levels Pw — 1W, no trasmis- 
sion is observed and the beam is rapidly absorbed along 
propagation. No significant time dynamics are observed. 
For co = 2%, at gelification, for power levels of lOOmW, 
a propagating light filament is observed, tunnelling on 
a time scale of the order of tens of seconds (Fig. [5]). 
The fluorescence signal is shown in fig. [5j The process 
is well rationalized by the theory and numerical experi- 
ments above. The relaxation time after the shock phe- 
nomena is in general longer than that measured by dy- 
namical light scattering or z-scan LJ], which is of the 
order of few seconds. This due to large structural mod- 



ification (e.g., enhanced packing) induced by the shock. 

Conclusions — We reported on a novel regime of 
optical-propagation in soft-matter, mediated by many 
effects: diffraction, absorption, temperature driven de- 
focusing and thermophoresis. The excitation of an hy- 
drodynamic shock, resulting from a dynamically arrested 
phase, allows propagation inside an highly absorptive sys- 
tem. This light-induced cavitation appears as a general 
process, demonstrating how complex processes underlie 
nonlinear wave propagation in self-organized matter and 
lead to novel spatio-temporal effects. 

We acknowledge support from CINECA. The re- 
search leading to these results has received funding from 
the European Research Council under the European 
Community's Seventh Framework Program (FP7/2007- 
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FIG. 5: (Color online) Experimentally retrieved fluorescence 
signal due to the formation of a propagating beam in the 
material at three different instants. The inset shows a three- 
dimensional representation at t = 100s (arbitrary scale). 
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